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Abstract 

We study the dynamics of Bose-Einstein condensed atoms in a 1-D optical lattice potential 
in a regime where the collective (Josephson) tunneling energy is comparable with the on-site 
interaction energy, and the number of particles per lattice site is mesoscopically large. By directly 
imaging the motion of atoms in the lattice, we observe an abrupt suppression of atom transport 
through the array for a critical ratio of these energies, consistent with quantum fluctuation induced 
localization. Directly below the onset of localization, the frequency of the observed superfluid 
transport can be explained by a phonon excitation but deviates substantially from that predicted 
by the hydrodynamic/Gross-Pitaevskii equations. 
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Coherent control of the collective dynamics of macroscopic quantum systems is currently 
of great interest due to possible applications in quantum measurement and information sci- 
ence jl, 2]. For example, coherent manipulation of superconducting currents in Josephson 
junction circuits has led to the realization of high-Q electronic qubits for quantum logic de- 
vices [3]. Similarly, manipulation of the superfluid properties in atomic systems, such as with 
BECs in optical lattices, may soon provide a realization of de Broglie wave interferometers 
which perform below the shot-noise limit {4]. 

At zero temperature, the physical characteristics of coupled superconducting/superfluid 
reservoirs are determined by two competing energies: the kinetic energy associated with 
tunneling between sites (Ej), and the on-site interaction energy (Ec), resulting from (re- 
pulsive) inter-particle interactions . Specifically, for a BEC in an optical lattice system 
Ej = N'j, where N is the number of atoms in a lattice site and 7 is the inter-site tunneling 
energy, and Ec = gf3, where g(3 is the mean-field energy. The nature of the many-body 
ground state in the lattice array is governed by the ratio F = Ec/Ej = gf3/N'y, which 
can be divided into three regimes. For F ^ 1, the system exhibits global superfluidity and 
long-range phase order. As F approaches 1, interactions lead to a frustration of long-range 
phase order and a corresponding reduction of the single reservoir number variance 6N, with 
SN ~ (N'j/ gPy^^. When F ~ 1, the system undergoes a transition to an insulating regime 
where number fluctuations are strongly suppressed {6N < 1) for commensurate filling. For 
translationally invariant lattice arrays, this defines the Mott-insulating regime (MI) [7, ^. 
This system has been shown to map onto 1-D superconducting chains, which demonstrate 
a Kosterlitz-Thouless transition 91]. 

Previously, interferometric techniques have been used to study the ground state properties 
of BECs in optical lattices forming 1-D arrays with high filling factor in the regime F > 1 
[lo| . This work demonstrated that, for the system studied, interferometric measurements 
do not have the specificity to reveal possible abrupt changes in the many-body state of the 
system as the Mott-insulating regime is reached HI]. However, recent theoretical work has 
shown that transport measurements in the regime F ~ 1 are expected to provide additional 



insights into the superfluid properties of the array [12[]. Similar transport measurements have 
previously been used to study the superfluid properties of arrays in the semiclassical regime, 
F <^ 1, and have observed both superfluid transport jisl] as well as damping. Dissipation 



has been attributed primarily to either dynamical instabilities [l^ or quantum fluctuation 



2 



effects [l5|,Q, 

depending on the experimental parameters. Both of these mechanisms make 
it difficult to observe finite amplitude transport for systems in their ground state at deep 
lattice depths, even for F < 1. 

In this work, however, we observe residual coherent transport in this dissipative regime, 
which we attribute to the presence of anomalous phonon excitations. We study this trans- 
port through direct observation of atom motion across a 1-D lattice array superimposed 
on a weaker harmonic potential. We observe a crossover from coherent oscillatory behav- 
ior to localization as we vary the strength of the coupling between adjacent lattice sites. 
We attribute this localization to the role of quantum fluctuations in driving the formation 
of strongly correlated, mesoscopic insulating states with lattice site atom number occupa- 
tion, ~ 100. We support this interpretation by observing the cessation of macroscopic 
tunnelling to occur when F ~ 1, the same critical relation found in analogous Josephson 
junction experiments 



17j. Furthermore, we rule out localization due to known semiclassical 
effects. 

The Bose-Hubbard Hamiltonian is expected to provide an accurate theoretical description 



of a lattice array of bosonic atoms 181]. Written in terms of the on-site single particle creation 
and annihilation operators aj and Oj: 

H = -lYl + Yl ]:gPia\a\aiai + ^ Via\ai (1) 

<i,j> i i 

where Vi = Qi"^ is the external potential due to the harmonic magnetic trap, the subscript i 
denotes the i th lattice site, and g = Anh^a/m (a is the repulsive s-wave scattering length 
and m is the atomic mass). /3j and 7 are determined from integrals over single particle 
wavefunctions [lol, and can be precisely, experimentally controlled by raising and lowering 
the intensity of the optical lattice in order to vary F. 

In our experiments, we investigate the response of the lattice ground state to sudden 
shifts in the position of the harmonic potential used to initially create the condensate. 
Exact theoretical predictions using Eq. 1 for the ensuing array dynamics are difficult for 
our experimental conditions due to the large Hilbert space needed to model the system. 
Nevertheless, we gain intuition into the dynamic behavior by considering limiting regimes. 
For static potentials and translationally invariant arrays, a second order quantum phase 
transition from a superfiuid to a MI phase is predicted to occur at F ~ 1 



211. 
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At this point the many-body wavefunction localizes to a product of states with nearly 



quantized atom number at each lattice site for commensurate lattice filling. Such states lack 
macroscopic phases, and thus, we expect suppression of Josephson-like transport of atoms 
through the array when an external driving force is applied. On the other hand, just below 
the transition, particle fluctuations 1 < 5N <^ N^^'^ are sufficient to to define an average 
macroscopic phase and induce superfluidity, even for large filling factors. 

Inclusion of the external harmonic potential (as in Eq. 1) complicates the analysis 18|, 
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24| . For deep lattices, where F > 1 across the entire array, mean-field analysis 



predicts isolated Mott domains of incommensurate filling with local incompressible regions 
but does not demonstrate the vanishing global compressibility indicative of a true phase 



transition 



2J]. These domains have been recently observed in Refs. 25|, |26(]. For shallower 



lattices, regions near the edges of the array can locally fulfill the Mott-insulator condition, 
while the central lattice sites remain in the superfiuid regime. We expect the formation of 
these insulating barriers to substantially alter the dynamics of the superfiuid confined in the 
interior regions. 

For a given lattice strength, U, we determine the values of gjSi and Ni using the follo wing 
two-step approach. First, we numerically solve the 3-D Gross-Pitaevskii equation (GPE) 27| 
for a single lattice site (neglecting tunneling between adjacent lattice sites) to determine the 
spatial wavefunction associated with the individual lattice sites as a function of the strength 



of the lattice potential and the number of atoms in the we^ 



1. This allows determination 



of the chemical potential fii and effective value of gf3i |28l . 1291]. Analytic expressions for 
these quantities can be obtained in the limit where the kinetic energy is negligible. For 
our parameters, the analytic expressions agree to within 10% with numerical estimates that 
do not neglect kinetic energy. Next, we estimate the equilibrium distribution of the mean 
occupancy of each lattice site across the array by equating the chemical potentials associated 
with each site, subject to the constraint that the total number of atoms sums to a fixed value. 
We obtain analytic expressions if we neglect kinetic energy associated with tunneling, which 
is an excellent approximation for the lattice depths explored in this work. For U = 50 Er, 
the central lattice site contains roughly 120 atoms, with a 1/e radius of 19 lattice sites, 
where gjS/h ~ 17 rad/s and 'j/hr^ 0.6 rad/s. 

Accurate determination of F (Fig. [1]) further hinges on measurement calibration of atom 
number and lattice depth. We infer the total number of atoms in the trapped conden- 
sate from an absorptive image of the atomic array. We measure the lattice depth in three 
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independent ways. First, we determine the lattice depth through direct observation of the 
harmonic oscillation frequency in each well: we measure atoms lost from the lattice via para- 
metric heating as a function of the modulation frequency of a small amplitude perturbation 
to the position of the potential minimum. Second, we measure the period of Kapitza-Dirac 
diffraction by suddenly turning on the lattice [30]]. Finally, we directly measure 7 at weak 
lattice depths {U < 10 Er) by observing the amplitude of the oscillation induced by an ap- 



plied transient phase gradient 



3l| . In particular, a phase shift (quasi- momentum) 6(f) = 7i/2, 



induces an oscillation with an amplitude (in lattice units) equal to y 27/^2 [32j. These ex- 
perimental techniques determine the lattice calibration with an uncertainty of ±8%, which 
agrees with the calculated lattice depth based on measured parameters to within 20%. 
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FIG. 1: 

r = Ec/Ej = g(3/N^ is calculated as a function of lattice depth. For a translationally invariant 
lattice system, F ~ 1 denotes the onset of a Mott-Insulating regime. 



Our experimental apparatus is described in detail in Ref. 101] and begins with 10^ laser 
cooled ^^Rb atoms loaded into a time-orbiting potential (TOP) trap. We use both TOP 
and forced radio frequency evaporative cooling to generate nearly pure condensates with 
4±1 X 10^ atoms, optically pumped into the F = 2,mp = 2 state. We estimate T/T^ ~ 0.27 
before the lattice is turned on (T is the temperature and Tc the BEG transition temperature). 
After the BEG phase transition, we adiabatically relax the magnetic trap, thereby increasing 
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the condensate's size and decreasing its density. This results in a radial trapping frequency 
of uj_[_ = 70 rad/s and Uz = 195 rad/s, corresponding to Q/h = 4.2 rad/s {Q = 2.1 x 
10~^i?/j). A retro-refiected laser beam in the vertical direction at a wavelength A = 840 
nm is superimposed over the condensed atom cloud to create a 1-D lattice potential. The 
light is focussed to a 1/e intensity radius of 50 fim and is far-detuned to the red from the 
780 nm optical resonance of ^^Rb. This creates a standing wave which forms a periodic 
array of potential wells spaced by A/2. Atoms occupy lattice sites at the anti- nodes of 
the sinusoidally varying optical potential. The lattice also provides additional transverse 
confinement of the atoms. For example, at f/ = 50 Er, the characteristic trap frequencies 
are 842 rad/s (transverse) and 3.1 x 10^ rad/s (longitudinal). 

We slowly increase the intensity of the lattice laser (in a time of 200 ms) to minimize non- 



adiabaticity in the state preparation. We note that 
for the deeper lattice depths used in this work 



;his ramp may not be strictly adiabatic 
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34| . In order to study the transport 



properties of atoms through the array, we suddenly shift the harmonic potential minimum 
and stroboscopically observe the ensuing motion of the center-of-mass of the array (using 
destructive absorptive imaging techniques). We shift the harmonic potential by pulsing on 
a weak vertical magnetic field. The resulting translation. A, to the minimum induces a 
chemical potential offset, E = AAQ/X, between adjacent wells, which drives subsequent 
array dynamics. 

We investigate three regimes of lattice transport: F ^ 1, F < 1 and F > 1. In order to tie 

rfn 

in with previous work p^, 115[, we first explore the superfiuid regime for low lattice depths 
(F ^ 1) where the semiclassical hydrodynamic equations are well suited. Transport oscil- 



lation frequencies closely follow a scaled harmonic magnetic trap frequency, oo = ujQ^Jm/m* 
with an effective mass, m/m* = 77r^/2i?/j, determined by the lattice depth 35|. We charac- 
terize the observed oscillations through their amplitude and frequency, determined by non- 
linear least squares fits to the oscillation time sequences [z{t) = ^e"*^/^™* cos(co't), where A 
denotes the oscillation amplitude and b the damping coefficient]. In Fig. [2]A, we see good 
correlation of our observed frequencies with semiclassical theory. Ground state coherent 
transport in this regime is also contingent on the trap displacement being sufficiently small 
so that an unstable regime is avoided. 

If the trap displacement is too large, the system enters a dynamically unstable regime, 
where small perturbations around a plane wave grow exponentially in time [l^. This is 
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FIG. 2: 

(A) Center-of-mass oscillation frequency vs. lattice depth for A = 3.5 ^m. Solid line is solution to 
semiclassical hydrodynamic equations. (B) Damping of center-of-mass oscillation vs. lattice depth. 



6/60 notation is adopted for comparison with Ref. 
harmonic trap frequency. 



151 ]. where 60 = 2mujQ, and loq is the magnetic 



expected to occur after a critical displacement, A„\t = X/2^j2'y /il, where a 7r/2 phase 
difference accumulates between adjacent wells [32]. In this regime, dissipation is expected 
to lead to rapid frustration of coherent tunneling. For A = 6.4 /xm, the largest displacement 
used in this work, this occurs for U > 10 Er (see Fig. [3]). 

The data shown in Fig. [21 however, demonstrates significant damping, despite being taken 
with A < Acrit- Furthermore, for this data 7 3> -E, indicating that the width of the lowest 
band is greater than the energy offset induced between wells (see Fig. [3]). Thus, we also 
do not expect localization related to Bloch oscillations. One possible explanation for the 
damping seen in Fig. EB is effects due to quantum fluctuations, predicted to occur even at 



low lattice depths [23|, 



361 ]. This effect has recently been suggested as a potential mechanism 
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FIG. 3: 

Critical displacement, Acrit, as predicted by the GPE, is plotted as a function of lattice depth, 
shown with a solid line. Above the solid line this semiclassical picture predicts a dynamically 
unstable, dissipative regime, and below the line, superfluid transport (SF) is expected. The dashed 
line indicates the condition where ^ = E. For E > ^ and weak interactions, lattice dynamics are 
governed by Bloch oscillations. 



for experimentally observed damping in Ref. [15!]. We note that the damping observed in 
Ref. [l5| is significantly greater than what we observe, however, in that system, lattice site 
occupation is on the order of unity. Other possible explanations for our observed damping 
include finite temperature effects [37 1. 

For r < 1 we observe higher frequency oscillations which emerge on top of the overdamped 
slower oscillations. Fig. IHA. shows a slow oscillation in the underdamped region with uj = 
88.2 ± 1.9 rad/s at f/ = 7.2 En. Fig. HfB displays the emergence of a faster oscillation with 
UJ = 149.9 ± 5.8 rad/s, taken at U = 10.8 Eji, where the non-linear fit requires an added 
linear term to account for the initial slope. We continue to observe high frequency transport 
oscillations as the lattice is increased for U < 50 Er, as shown in Fig. [5JA.. We rule out the 
possibility of these high frequency oscillations being Bloch oscillations since u is independent 
of A, which determines the energy offset between adjacent wells. 

The presence of these oscillations cannot be explained within the ground state semi- 
classical picture, since the observed frequencies are much greater than those expected from 
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FIG. 4: 

Center-of-mass oscillations with A = 3.5 //m are shown for (A) the underdamped regime, with 
uj = 88.2 lb 1.9 rad/s for U = 7.2 Eji and (B) the overdamped regime, with uj = 149.9 it 5.8 rad/s, 
taken at U = 10.8 Er. 



the hydro dynamic equations (see dashed hne in Fig. 04.) [35|. In addition, our observed 
frequencies cannot be explained by using standard approximations to include the role of 
quantum fluctuations in ground state transport. Recent theoretical approaches using the 
Truncated Wigner approximation (TWA) have shown that quantum fluctuations tend to 
reduce transport frequencies in the overdamped regime and to increase damping 2^. 

We expect, however, that the observed frequencies in this regime should be related 
to the characteristic frequencies for phonon excitations. An array phonon frequency can 
be determined for an effective phonon wave- vector, gefr, which characterizes the array co- 
herence length. We infer geff from the phonon dispersion relation for a uniform lattice 
hug = ^Ji-f sm'^{^)[2Nigl3 + A-f sin^i^)] [d]. We find that for suitably deep lattice depths 
(near the expected MI cross-over), our observed frequencies correspond to phonon excita- 
tions with 47r/geffA = 4 lattice sites, as seen in Fig. [5jA.- It is interesting to note that this 
phonon frequency also corresponds to the generalized Josephson frequency associated with a 
2-well system. As described earlier, the complexity of the Hilbert space of our system makes 
achieving a quantitative theoretical prediction for ggfr very difficult. However, our empirical 



observation of an effective phonon excitation length (4 sites) , much less than the full spatial 
extent of the atomic cloud (19 sites), may be due to the emergence of MOT domains in the 



outer wells 
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261]. For shallower lattice depths (in the regime shown in Fig. 0\) we observe 



frequencies which are consistent with longer range phonon modes, possibly due the effective 
array coherence length increasing with reduced F [38]. 

We note that for these lattice depths, our experiments probe the ill-understood regime 
where the semiclassical (Gross-Pitaevskii/hydrodynamic) equations become dynamically un- 
stable (see Fig. [3]). In this regime, for example, the Bogoliubov dispersion relations have 
imaginary solutions and unconventional normalization 39|], and thus, the TWA is unreli- 
able. It has been shown that the exact quantum dynamics of a dynamically unstable system 



diverge (logarithmically with A^) from the classical trajectories 40|, |4l|, which further com- 
plicates the analysis of our excited state transport. 

For F ~ 1 we observe an abrupt cessation of superfluid transport. The oscillation ampli- 
tude falls to zero at a critical value for the lattice depth corresponding to F ~ 1 (Fig. [SB), 
consistent with the expected superfluid-MI transition [ol, [21I (at f/ = 49 E^, F = 0.21 with 
an error range of 0.09 < F < 0.52 determined, predominantly, by our experimental uncer- 
tainty in lattice calibration). Fig. [5]B (see lower inset) illustrates the oscillation induced by 
a shift (A = 6.4 fim) in the harmonic potential at a lattice depth of U = 47 En and is con- 
trasted with localization observed for the same displacement at U = 54 Eji (see upper inset). 
The observed threshold is independent of the lattice displacement (within our experimental 
limits). Additionally, above the transition, the atoms are observed to localize immediately 
after the shift in position of the harmonic potential, as illustrated in the inset to Fig. [5B- 



These observations exp 
quantum self-trapping 



icitly rule out other localization mechanisms - such as macroscopic 



42, 



43| or dynamical instabilities 32| - where the localization de- 



pends on the initial displacement, or manifests itself only after the array dynamically evolves 
following the displacement (in the case of dynamical instabilities). 

In conclusion, we have used the presence of anomalous phonon excitations to probe 
atom transport deep in the dissipative regime. This regime has been previously inaccessible 
using ground state array transport. We observe localization to occur at the critical value 
F ~ 1, which suggests localization due to quantum fluctuation effects. This work raises 
the possibility for further control of macroscopic quantum coherence, even in excited state 
systems. 
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FIG. 5: 

(A) Center-of-mass oscillation frequency, uj, vs. lattice depth for harmonic trap displacements A = 
6.4 /im (circles), 4.3 /im (triangles), and 3.2 /im (squares). The solid line indicates the excitation 
frequencies for a lattice array phonon mode with an effective length of 4 sites, and the dashed line 
denotes the oscillation frequency predicted by the semiclassical hydrodynamic equations. Filled 
points, shown with u = 0, represent data where we cannot experimentally resolve an oscillation. 

(B) Center-of-mass oscillation amplitude vs. lattice depth. Typical oscillation sequences are shown 
in the upper and lower insets, for parameters U = 54 En, A = 6.4 ^m and U = 47 Er, A = 6.4 
/im respectively. Below the transition, errors are inferred from the residuals of the non-linear least 
squares fits. Above the transition, error bars are determined from the variance of the ensemble of 
data points. Amplitudes are scaled to the measured oscillation amplitude for U = 42 Er. 
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